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Abstract. The methods of effective field theory are used to study generic theories of 
inflation with a single inflaton field and to perform a general analysis of the associated 
non-Gaussianities. We investigate the amplitudes and shapes of the various generic 
three-point correlators, the bispectra, which may be generated by different classes of 
single-field inflationary models. Besides the well-known results for the DBI-like models 
and the ghost inflationary theories, we point out that curvature-related interactions 
may give rise to large non-Gaussianities in the form of bispectra characterized by a 
flat shape which, quite interestingly, is independently produced by several interaction 
terms. In a subsequent work, we will perform a similar general analysis for the non- 
Gaussianities generated by the generic four-point correlator, the trispectrum. 
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1. Introduction 



The inflationary paradigm is central in modern cosmology as it naturally provides 
an explanation for long standing issues as the flatness and horizon problems. The 
simplest standard, single-field slow-roll models of inflation already nicely account for 
the scale invariant primordial power spectrum of the Cosmic Microwave Background 
(CMB) anisotropies [H El |3l IH O E]. On the other hand, in order to get a deeper 
understanding of the inflaton dynamics, one needs to consider observables sensitive to 
deviations from Gaussianity, which might be studied in the form of the three-point 
correlator of cosmological perturbations, the so-called bispectrum [TJ El [TOl [11], the 
connected four-point correlator, the so-called trispectrum [12], loop corrections to the 
power spectrum [13] and so on. In recent years, a lot of progress has been done in 
this direction resulting in the characterization of a wide variety of inflationary models 
through their higher order correlation functions [H]. 

Complementing this theoretical effort, new experiments provide a better sensitivity 
to deviations from Gaussian statistics (see, e.g. [171 [IE])- The very recent launch of the 
Planck satellite [T51IT6] (and the continued analysis of WMAP data [6]) presents us with 
the exciting opportunity to actually test this zoo of possibilities and makes more urgent 
the classification and refinement of the different predictions at the level of bispectrum 
and trispectrum of curvature perturbations (see, e.g., [191 [201 [21] for comprehensive and 
updated reviews). 
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In this context, the effective theory approach proposed in [22], [23] and further 
analyzed in [23] represents a powerful tool. Indeed, the effective Lagrangian approach 
provides an efficient way to generically study large classes of inffationary models. 
Furthermore, this formalism provides a clear-cut dictionary between deviations from the 
standard scenario (almost Gaussian perturbations) and higher-order invariant operators 
in the action; it sheds new light on effects due to symmetries in the action (notably, a 
reduced speed of sound often automatically results in an enhanced non-Gaussianity); 
it allows some general and immediate calculational advantages, by exploiting the so 
called decoupling regime; finally, one can classify almost any specific inflationary model 
switching on some particular operators in the Lagrangian, thus allowing a more unifying 
approach. The effective field theory formalism naturally comprises Lagrangians with 
higher order terms. These terms call for an UV completion of the underlying theory. In 
this perspective one looks with particular interest at string theory inspired models, such 
as DBI infaltion [27], and possible explanations in this context for k-inflation models 
[251126]. 

They are retrievable in the effective theory approach as the effective field theory 
that lives in a particular region of the parameters space. 

In the spirit of effective theory, one should also include all possible curvature- 
generated terms {5R^ ^,...) in writing the most generic effective theory Lagrangian of 
inflation up to the desired order in perturbation theory. In this paper we perform a 
thorough analysis of the amplitude and shapes of the three-point correlators which 
arise from the various classes of inflationary models driven by a single scalar field 
described within the generic effective theory approach. Besides reproducing the known 
results exisiting in the literature, that is to say that large non-Gaussianities may be 
generated if the sound speed is smaller than unity, we also show that single-field 
models of inflation may produce large non-Gaussianities in the form of three-point 
correlators whose amplitude is peaked in the so-called flat configuration, and precisely 
corresponding to a shape of the bispectrum where two of the momenta are roughly one 
half of the third one ki = 2/c2 = 2^3. Such bispectra emerge from curvature-generated 
interactions that have not been discussed so far in the literature. Let us recall that a 
shape of this type, peaked on flat triangles, was first found in [38] where it was obtained 
as the result of a linear combination of two operators with proper weights |. We stress 
here that our results further expand the space of operators that generates this kind of 
shape and allow one to easily obtain it from different, independent operators. Before 
the results of [38] , the only scenario compatible with such a flat shape (also called folded 
or squashed) for the primordial bispectrum was that of models of inflation which relax 
the assumption of a regular Bunch-Davies vacuum, as studied in [HI [281 129]. In fact, 
up to [38], it was standard procedure to associate the primordial bispectrum evaluated 
in the Bunch-Davies vacuum to the so-called equilateral or local shapes [10] •§ Here we 

I It is important to note that these weights are not fixed, indeed there is an order one interval in |38| 
that produces shapes compatible with the flat template. 

§ We refer to [30l [31] for a detailed analysis of possible shapes of the primordial bispectra and to the 
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highlight that such a flat configuration for the primordial non-Gaussianities might be 
generated in single-field models of inflation in a more general context than previously 
thought. It is worth pointing out another aspect in which our results differ from previous 
ones. We are able to provide a solution for the equations of motion of the primordial 
perturbations that can interpolate between known cases, like DBI or Ghost inflation. We 
are able to study in full generality these intermediate situations also when interaction 
terms are included. In a follow-up paper [32], we will perform a similar analysis for the 
trispectrum. 

The paper is organized as follows. In section 2 we give a brief account on how to 
write the Lagrangian we employ in the following sections, following Ref. [23]. The reader 
who is familiar with this approach might wish to skip this part and jump directly to the 
third-order Eq. ( fT3l) . We also comment on the inflationary models one can span within 
this approach and on the claim of its generality. The calculational algorithm used in 
this section can be extended to higher orders as well. In section 3 we derive the equation 
of motion for the scalar fleld which in the effective theory approach is used to describe 
the cosmological perturbations [33] and we provide an analytical solution to this general 
equation. In section 4 we present the calculation of the amplitude of the bispectrum 
for each interaction term in six different conflgurations capturing the simplest single 
fleld slow roll model, DBI inflation, ghost inflation and intermediate cases. We focus on 
regions in the parameter space where non-Gaussianity is large. In section 5 we study 
the shape of the bispectra arising from the individual interaction terms identifying, 
among other things, interesting shapes from the curvature-related contributions. We 
then summarize our flndings and comment on further work in section 6. 

2. The inflationary action up to third order in the eflPective theory approach 

Our initial goal is to write down the complete theory of single-fleld models of inflation 
up to third order in perturbations. We will follow the effective theory approach flrst 
introduced in Ref. [23]. Let us give a brief summary of the procedure. The scalar fleld 
responsible for inflation is splitted as usual as an unperturbed part plus a fluctuating 
one: 



One then chooses to work in the the comoving (or unitary) gauge for which 6(f) = 
[22] . The Lagrangian will no more be invariant under full spacetime diffeomorphisms 
(diffs) but only under the spatial reparametrizations. This is the starting point to 
write the most general space diffs invariant Lagrangian at the desidered order in 
perturbation theory in an effective theory approach. In [23] the authors prove that, 
once an approximate shift-simmetry is operated, their second and third order action is 
the most general one (see also [23] for an interesting perspective on the most general 
effective Lagrangian for inflation). One can then use the Stueckelberg procedure to 

discussion in Ref. [38] for an approach that includes the flat shape, (see also Ref. [TT]). 




(1) 
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restore full spacetime invariance. As a by-product of this procedure, the degree of 
freedom hidden in the metric shows up again as a scalar field. 

The most general space diffs-invariant action in unitary gauge is schematically 
written as 123): 



S 



9 F{.R^upcn 9 1 -^^ly^ 



(2) 



where i^^j/ is the extrinsic curvature on which we will elaborate more later and the 
indices on g^^ are free indices. Considering fluctuations around a FRW background 
amounts to studying the following action 



S 



rl 



-9 

„00 



n>2 



(3) 



where the -F*^"^ functions contain fluctuations which are at least quadratic. The next 
step is to restore full spacetime reparametrization invariance. To see how it works, we 
borrow a simple example from [23] and consider the following sample action terms 



d^x A{t) + B{t)g' 



00, 



x\ 



We are interested in time reparametrization t ^ t + ^oi^i t)] 
above action (after a simple variable redefinition) reads 



(4) 

X — X, under which the 



g^^{x) 



(5) 



dx^ dx'^ 

Upon promoting to a field, ^o{x) = —7i{x) and requiring the following gauge 
transformation rule 7r{x) — )■ 7r(x) — .^0(3^) on it, the above action is invariant under full 
spacetime diffeomorphisms. The scalar degree of freedom vr makes its appearance in the 
time dependence of the A, B coefficients and in the transformed metric. Furthermore, 
under time reparametrization the metric g^^ transforms as follows: 

dg°' dg^ 



9 



(9x' 



:9 



which implies 



9'' -^9''; 
9'' (1 



vr) V° 



(6) 



(7) 



'l + rr)g'' + g^^d,7T; 
^2{l + n)g'^%7i + g'^di7id,7i. 

This procedure has been borrowed, conceptually unchanged, from standard gauge 
theory: a Goldstone boson which transforms non linearly under the gauge 
transformation provides the longitudinal component of a massive gauge boson. At 
sufficiently high energy such Goldstone boson becomes the only relevant degree of 
freedom. This is the so-called equivalence theorem. The same is true for our case: 
for sufficiently high energy the mixing with gravity becomes irrelevant and the scalar 
vr becomes the only relevant mode in the dynamics. This is the so-called decoupling 
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regime. Let us clarify this concept with a simple example. Consider the following 
contribution, taken from Eq. ([3]) 

M^iif^oo ^ Ml^H{{l + 7r)2^°° + 2(1 + ii)g'''diTi + g'^ diTidjTi). (8) 

We focus on the quadratic part of the first term in the above equation. Upon canonical 
normalization, tt^ = Mpiif^/^Tr and g^ = Mpig'^'^, one gets 

M^,H{g^' + 27ig'^ + rr^g'^) = + 2H^lH,gl'^ + Hg^^. (9) 

Consider the second term of Eq. ([9]) which mixes gravity with the scalar. Since tTc ~ Euci 
at energies higher than ~ H^^^ the term tt^ dominates the dynamics. This turns out to 
be true in general: the number of derivatives (which in Fourier mode would basically 
give an energy-dependent coefficient in front of tt) is higher in terms containing only vr's 
than in the mixed terms and therefore there exists an energy threshold above which the 
scalar decouples from gravity. Since in explicitating the F^'^^ term in Eq. ([3]) there can 
be, in principle, other quadratic terms that go like tt^, one has to consider which one 
is the leading kinetic term and determine the canonically normalized field tTc and the 
energy threshold accordingly. To take the safe route, one might well take the energy 
threshold, -Emixj to be the highest one of this set. Since one is concerned with correlators 
just after horizon crossing, one concludes that the decoupling procedure works as long 
as the decoupling energy is smaller than the Hubble rate H . More precisely, we can 
anticipate that the kinetic terms in F*^^) which are going to matter in our discussion 
come with coefficients M| and Mp^eH"^. The condition -Emix < -f^ is then satisfied if 
Mp^eH^ > M|, where e = —H/H^ is a slow-roll parameter; if this is not the case we 
need to assume M| < Mp^H"^. 

From now on we will work in the decoupling regime. In considering the terms of 
Eq. ([3]), we will therefore use only the unperturbed entries of the metric tensor. To 
write the effective Lagrangian up to third order, we start from Eq. ([3]) and follow the 
algorithm given in [23]. Fluctuations are encoded in the F^""^ terms. In order to be as 
general as possible, we also include all possible contributions up to third order coming 
from extrinsic curvature K^^ terms. In fact, it is instructive at this stage to step back 
and consider the action in Eq. ([2]). Given a theory which is space diffs- invariant, one 
can always identify a slicing of spacetime, described by a timelike function t{x), which 
realizes time diffeomorphism: on surfaces of constant t the time symmetry breaking 
scalar is also constant. Before selecting a gauge, there is still the freedom to make a 
choice on t{x) and working in the unitary gauge amounts to requiring t = t. In order 
to describe the geometry of this preferred slicing, one employs the extrinsic curvature 
tensor. In writing down such a tensor, one needs two ingredients: the unit normal 
vector n^, perpendicular to the constant t surfaces and the induced metric h^y. These 
are defined as 



; = g^y + n^n^, (10) 
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which allows us to write 



2(_^00)3/2 ' 2(-^00)5/2 2(-^00)V2 ' ^^^"^ 

The above expressions can be used to write explicitly the most generic third order action 
for the fluctuations around the FRW background within the effective theory approach 
ofRef. 123]: 



-9 



'^M^,R + Ml,Hg 



00 



Ml,i3H' + H) + ^M,it)\g'' + If 



2 
3! 
3! 



5K^'JK%5K% 



Mr{t) 



3! ^ 



3! 
3! 



(12) 



The coefficients M2, M3 and Mi, ■ ■ ■ , Mg are to be considered generic. The M coefficients 
multiply curvature-generated interactions. A given particular set of them will specify a 
given inffationary theory. The action f|T2|) is not yet invariant under full diffeormophisms 
though. One needs to follow exactly the steps illustrated in Eqs (jl]), (l5j) and (IE]) and 
promote ^0 to a field vr with the proper gauge transformation. 
In the decoupling limit we find: 



d Xy/—g 



M,{t? f idh)id]7r) + 2H{d.7Tydfn + 27T^,^.7r^,7^ ^^ 2 1 . M.jtf tt 



+ 



a" 



a" 



,4 ^"'^J 



3! a6 3! a6 



Mg{t) 1 



3! a' 



^dijirdjkTrdkiir 



A few clarifying comments are in order: 



If we consider terms only up to second-order, for M2 = Mi,2,3 = one recovers the 
usual quadratic Lagrangian for the fiuctuations, with sound speed = 1 and the 
standard solution to the equations of motion. Switching on M2 amounts to allowing 
models with sound speed smaller than unity, l/c^ = 1 — 2M|/(Mpiif), which are 
often linked to a high level of primordial non-Gaussianity [HI |23]. Furthermore, 
turning on M2,3 in the de Sitter limit, one recovers Ghost infiation [35]. On the same 
lines, keeping all the M's vanishing, but going to third and higher order with the 
M's, one can retrieve the interactions that describe DBI infiation [271 El EZ] • The 
list of correspondences continues with K-infiation theories and others, thus showing 



(13) 
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how the effective action approach provides a unifying perspective on inffationary 
models [23] . 

• The action in Eq. (fT3|) has already been written with large non-Gaussianities in 
mind. This means that, at every order in fluctuations and for each M and M 
coefficients, we have selected those leading terms which will eventually generate 
large three-point correlators. To clarify this point, we provide a simple example. 
Let us consider the terms up to second order in Eq. ( |T3i) and set conveniently 
^1,2,3 = 0. The properly normalized solution to the equation of motion will be 
the usual 7rfc(r) oc iHe~^^'^"'^{l + ikcsr). It is straightforward to verify that, at 
the horizon crossing , tt ~ Hn and Vvr ~ H/cg tt. Therefore, among the it terms 
with the same number of derivatives, the ones with the highest number of space 
derivatives dominate in the ^ 1 limit. Generalizing these estimates for the 
classical solution (which we will describe below) obtained from the equation of 
motion of our complete action, one selects the terms in Eq. (fT3l) . 

• there is also the comparison between same perturbative order but different M terms 
to be made. In the literature, all non zero coefficients in front of the various 
operators are generically assumed to be of the same order (see for example the 
discussion concerning the orthogonal configuration in Ref. [3H] for an interesting 
perspective). We shall not restrict ourselves to this situation. Note that, were the 
coefficients to be all of the same order, one could already identify the dominant 
operators. For example, consider the third order contributions ((9j7r)^(9|7r/a^ ~ 
M^[H'^/c^)7T^] for ^ 1 this will be a leading contribution with respect to, 
say, M2{tfH{dlTT){djTif ~ H {H^ / cfjTx^ . This is due to the fact that in the 
effective Lagrangian every additional derivative comes with a H/M <^ 1 factor 
attached: one is basically doing an H/M expansion where M is roughly the energy 
range of the underlying theory. In the last example we have intentionally picked 
terms with the same power of Cg at the denominator. Let us now look at the 
M'^{t){df'n')^ ~ M{H^ /c^)tt^ term though; comparing this contribution with the 
Ml term amounts to comparing M^ with MH'^/cl. We see that for a very small 
speed of sound the M7 contribution may still be relevant. These examples justify 
our strategy of including all the terms in Eq. (fT3!) compatible with ^ 1 (we will 
make more comments on this point in Sec. HI). 

Having written the action, we now proceed to elaborate about the solution to the 
equation of motion for the scalar degree of freedom tt. 

3. Solution to the equation of motion for the scalar degree of freedom 

To solve for the equation of motion for the scalar degree of freedom tt, we write the 
Lagrangian to second order in the perturbations. Note that in the equation below, as 
opposed to what was done in writing the third order action in Eq. (fT3l) . we include all 
possible terms, even the ones which are subleading and are therefore not expected to 
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give large non-Gaussianities. 



£2 = a 



3i/^^ + ^(5^)^ - ^ 3i^V - if^^ + -AdM' . (14) 




Following a common procedure, we have suppressed the time dependence in the M 
coefficients since we choose to work at leading order in the generalized slow roll 
approximation. After the usual change of variable, 7r(fc,t(r)) = u{k,T)/a{T), the 
equation of motion reads 
2 

u" -u + aok'^u + Pok^r'^u = 0, (15) 

where a^, f3o are time independent (again, at leading order in slow roll) dimensionless 
coefficients: 

-M|iij - MfH/2 - 3/2 M^H^ - M^H'^ /2 



-Ml^H + 2M| - 3Mf /J - 9/2 M^H"^ - 3/2 MlH"^ 



n ^ (Ml + Ml)H' ^g^ 

2{-M^^H + 2M| - 3MfH - 9/2 M|iJ2 .3/2 Mfi/^) ' ^ ^ 

Now, for the sake of simplification and in order to make sharper and more clear 
the correspondence between switching on the coefficients and recovering any specific 
inflationary models, we make the following choice: we take M| = — 3M| and introduce 
Mq = 2/3 M|. As a consequence, the second order Lagrangian now reads: 



= [M^Md.Trr + 2M^e - M^H ^3^ - - ^ [-.i^t^f)) ■ (17) 

Since in the Lagrangian we lost one coefficient in going from M2,3 to Mq one might well 
count one less degree of freedom at our disposal in the calculations that follow. On 
the other hand, as one can easily check, the operators with the M2, M3 coefficients give 
rise to almost the same amplitude and shape for non-Gaussianities. Thus, the price to 
pay is not that high. The ao and /3o parameters generalize the sound speed and the 
H'^/M'^ term (ffist introduced in Ref. [35]), respectively. Indeed, considering from now 
on Eq. (ITTj) as our starting point, the complete expression for these parameters reads: 

-Mgiij - MfH/2 MlH' 



-M^^H + 2Mi -3MfH 2{-M^^H + 2Mi - 3MfH) 

so that, for instance, one recovers the sound speed 1/c^ = 1 — 2M|/(M|[if) for Mi = 0. 
The solution for standard inflation is recovered if Pq = Mi = 0; the one for ghost 
inflation is obtained in the de Sitter limit if ao = = Mi. Notice that, as evident 
from Eq. f|T8l) . there is the possibility of a negative value for oq, the generalized speed 
of sound. By looking at Eq. (|T5|) one sees this would immediately result in a dispersion 
relation which, for aok"^ dominating over /3ofc^r^, gives an exponential solution for the 
modes. As explained in [38], this fact does not necessarily exclude models with negative 
ao. In fact, as far as the /3ofc^r^ term dominates soon enough in the UV before the 
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cutoff scale of the theory is reached, these models are still in principle acceptable as 
they eventually freeze at the Hubble scale. On the other hand, having an exponential 
growth, even for a limited energy window, can generate large non-Gaussianities which 
are ruled out by observations [38]. We plan to perform a detailed analysis of the issues 
associated with the ao-negative region of the parameters space in [39] . 

We can readily make some educated guesses on the behaviour of the wavefunction 
before solving the equation of motion. First of all, the typical oscillatory behaviour deep 
inside the horizon is to be expected: in this regime both aok"^ and fSok^r^ cause wave-like 
behaviour on the wavefunction while the (— contribution is negligible. This tells 
us the main contribution to correlation functions will be coming, as usual, from the 
horizon crossing region and the outside. Note here that, if Pq ^ 0, the 'ghost inflation' 
term will eventually lead the oscillation if one goes deep enough inside the horizon. On 
the other hand, in the r — )■ limit, the term leads the dynamics and we recover the 
usual frozen modes. As done for the DBl case, it is convenient to introduce the notion 
of an effective horizon crossing, placing it where the oscillatory behaviour stops being 
dominant 

aok^ + (3ok'T^ ^ r, = ^ (19) 

^* ^Y"o + \/ao + 8/^0 

For /3o = and ao ~ 1 one recovers the usual k'^T^ ~ 1 horizon crossing. 

Let us now be more quantitative. The solution to Eq. (fT5|) . being of second order, 
will come with two momentum-dependent integration constants. We have determined 
their values by requiring to obtain the known standard and ghost solutions in the 
corresponding limits. The general wavefunction reads 



u{k, t) 



21/4, 



-26 2 



HyperG 



1" Vl'"2' 



iJ/3ok'T' Ci{k) 



+ LagL (\ + -| -zV^A^VM C,{k) 



(20) 



where HyperG is the confluent hypergeometric function and LagL stands for the 
generalized Laguerre polynomial. Upon fixing 



Ci{k,aQ,f3o) 





ao + VA)^ 


-3/4 

r 


5 iag 


J.-3/2 




- 3M^H T 


3/2 

4 







C2{k) = 0, 



one can recover the normalized wavefunctions of standard inflation [36] and ghost 
inflation [35]. We can rewrite the solution as 



i{k 





5 iao 


HyperG 


[ ^ 


4^' 2' Wl^okrj 


^2(M|iei/2 + 2M| - 3MfH) (ao + 




'"k^nv 


5 1 ao 
_4 4ao-4ii/AD. 



(21) 
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whose asymptotic expression for r — )> reads 
Tr(k, r ~ 0) ~ 



Ml^em + 2M| - 3Mfi7 (a^ + ^of'^ k'^/^ T 



(22) 



Note here the normahzation factor in the denominator of Eq. (|22|) : 



'Ml^em + 2M| - mfH = ^^. 

By looking at Eq. ( IT3|) one understands that this factor originates from canonically 
normahzing the kinetic term, which leads to to tt^ = A/'tt . Let us also mention that 
the expected behaviour of the wavefunction at horizon crossing we commented above, 
77 ~ Hit and Vtt ~ H/cg vr, are now generalized to 

V2H 



TT ~ 



Hn, 



Vtt 



TT 



(23) 



"0 + 



an 



It is then clear how one can in principle obtain high level of non-Gaussianity by means 
of terms in the Lagrangian with space derivatives (requiring ao, /3o 1 in the process) 
rather than time-like ones. A more detailed analysis of the general solution to the 
equation of motion and of the resulting power spectrum is to be found in [39] . 



4. Amplitude of the primordial non-Gaussianity 

In this section we wish to perform a general analysis of the amplitude of the bispectra 
stemming from the general third-order interaction terms. The shape analysis will be 
done in the following section. In the calculations that follow we employ the so-called 
in-in formalism [iQl HU SSj |13]. To compute the amplitude of the non-Gaussianity, 
indicated by /nl, we proceed as traditionally done in the literature and evaluate the 
three-point correlator in the so-called equilateral configuration where all momenta are 
taken to be equal: ki = k2 = k^. In other words one can write the bispectra of the 
gauge-invariant curvature perturbation ( generated by each interaction term (I) as 

(C(ki)C(k2)C(k3))/ = (27r)353(ki + k2 + ks)Bi{k,, k2, ks) , (24) 

with an amplitude f^-^ defined so that with all three-momenta equal /j^l = 
{Q/5)Bj{k,k,k)/P^{k)'^, where P({k) is the power spectrum of the curvature 
perturbation. For large non-Gaussianities the linear relation ( = —Hit suffices for 
the bispectrum calculations since quadratic corrections give a negligible contributions. 
Notice also that the values of the integrals appearing in the in-in computations have been 
specified at horizon crossing and the contribution of the integral function at r = — oo 
has been put to zero mimicking the effect of the slight rotation of the r-axis into the 
imaginary plane. 

Given the broad number of possibilities, we choose to compute numerically the 
amplitude of the various bispectra identifying six benchmark points in the (ao, /3o) plane 
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and numerically integrating the exact wavef unctions. Since one is interested in probing 
models with large non-Gaussianity, the values of [oq, Pq) are taken much smaller than 
unity as described in the following Table: 



Benchmarks 


1 


2 


3 


4 


5 


6 


«o 


10-2 





0.5 ■ 10-2 


2 ■ 10-^ 


10-4 


10-6 


/3o 





0.5-10-^ 


0.25 • 10-^ 


5 • 10-5 









Upon using Eq. ( ITSll one can check that the first four benchmark points correspond to 
the same choice of the effective horizon. This choice has been made to suitably perform 
a comparison between the various cases and in particular against the values of /nl of 
purely DBI and Ghost models which correspond to benchmarks 1 and 2. The last two 
configurations (which are pure DBI) probe the space of extremely small ao and, in 
interactions with at least two space derivatives, are expected to give a larger amplitude 
than the first four points, at least if the interaction terms are regulated by unconstrained 
masses. To get the feeling for the figures involved, if we restrict ourselves to the case 
of theories for which ao reduces to the usual sound speed c^, typical values of are 
between 10-'^ and 10-^. Note also that the definition of Cg varies according to which 
operators are switched on in the action (one can well be in the de Sitter limit where 
the only spatial quadratic term has four derivatives; this leads to a different Cg, see also 
[38]). Therefore we choose to specify all the amplitudes as a function of the various M 
and M masses. 

It is convenient at this stage to clearly point out which M coefficients in Eq. (fT3ll 
are free and underline the relations among the constrained ones. From Eq. f|T8l) one can 
see that, despite fixing ao, as we did in the Table, there is still the freedom to pick any 
reasonable value for either M| or Mf. Similarly, fixing /3o does not completely specify 
Mq. In other words, both Mf and Mq are constrained by our choice of the (q;o,/3o) 
parameters; all the other coefficients are unconstrained. As the M's are expected to 
set the energy scale of the various underlying theories, they should be larger than the 
Hubble rate H, and can go up to Mpi. As elucidated in Ref. [21] though , for the action 
(fT3|) to be as general as possible, one might want to require the M's to be smaller than 
the Planck mass. That said, some useful inequalities that the mass coefficients must 
respect can now be reminded. Due to the fact we are working in the decoupling regime, 
we must require M| < Mp^H"^. Also, the fact that we are probing the (ao, /3o) ^ 1 space, 
imposes bounds on some masses. Consider the parameter ao in Eq. ( fTSi) . There are two 
ways this coefficient can be much smaller than unity. The first and perhaps most natural 
way, is to ask M| ^ Max {—MpiH, —MfH) which, due to decoupling inequalities on M2 
puts a bound on Mi, Mf -C Mp^H. The other possibility requires a partial cancellation 
in the numerator of ao, —Mp^H ~ MfH/2 which is certainly possible but it implies we 
are neither in the DBI, nor in the ghost regime, both of which have Mi = 0. Looking 
at (3q we see it is enough to require M^H'^ <^ M| or M^H"^ <^ —MfH and again, the 
first condition seems more natural. Let us stress here that, upon requiring the masses 
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to be all of the same order, M, and using that H/M <^ 1 in the effective theory, one can 
easily obtain small ao,/3o coefficient. However, when employing a single mass scale M 
in the whole Lagrangian, working with tiny values for and /3o would put a bound on 
M and necessarily influence the magnitude of all the interaction terms. In our analysis 
we let the M, M's coefficients be not all of the same order (with some important caveats 
upon which we expand at the end of this section). 

Below we present the results for each interaction term. All the amplitudes can 
be written as a dimensionless coefficient, 7„, times an (ao, /3o)- dependent numerical 
coefficient. The terms described in the first subsection are interactions that have already 
been discussed in the literature. The novelty here is represented by the fact we are able to 
study also interpolating configurations through the third and fourth benchmark points. 
In the second subsection we report on the amplitudes of the contributions from some 
curvature-generated terms that have never been discussed in the literature. 

4-1. Amplitudes from DBI-like interactions and first two curvature- generated terms 

The amplitudes from DBI-like interactions and first two curvature-generated terms are 
the following: 



-2M|7r(9,7r)Va^ 
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where 



71 



Ml 



2M| + Ml^eH"^ - ?,MlH' 



(25) 



We see the parameter 71 can in principle be of order unity. Indeed, if one assumes M| 
is the largest term in the denominator (in DBI this would correspond to a very small 
speed of sound), 71 is roughly 1/2. 



O2 = -4/3 MItt^ 
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where 



Ml 



2M| + M^^eH^ - 3MfH' ^^^^ 
The parameter 72 can be even larger than unity if, for instance, M3 is larger than M2. 
We can see, already at this stage, the effect of small values of ao and /3o at work: the 
numerical factor of a spatial derivative-free interaction is much smaller than that of a 
third order term like the M2 one calculated above which has two spatial derivatives. 
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. = -1/2 Ml{dinfd]n/a^ 
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Due to the fact that Mf is a constrained parameter we find that the amphtude equals 
an (ao, /3o)- dependent number times the paramter 

^ -Ml,eH\l-a,) + 2aoMl 

M^yeH^ + 2Mi -3MfH ' ^ ^ 

One should substitute the various values of ao in the Table in the 73 expression given 
above. Barring cancellation between different mass terms, 73 is generally smaller than 
one and can be as small as ao itself. 

• C»4 = -1/2 MoV4 {5H{df7r){d,7ry + ATidfd.ndjTr) /a^ 
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with 

74 = {MlH'^)/{Ml^eH'^ + 2M^ - 3M^H). (28) 
We see that large numerical coefficients appear in such a case, nevertheless 74 = /3o ^ 1. 

4-2. Amplitudes from curvature-generated novel interaction terms 

We come to the curvature-generated interaction terms that generate novel bispectra. 
Their amplitudes are given by 

• O5 = -2/3 Mi-jr^dfTT /a^ 
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where 

75 = {M^H)/{M^^eH^ + 2M^ - 3M'^H). (29) 

The coefficient 75 can be larger than unity. Barring cancellations in the denominator, a 
75 ^ 1 larger than unity imposes M4 ^ M2, Mi. 
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. = 1/3 Ml7r{dfn)^ /a^ 
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where 

76 = {M^H^)/{Ml,^eH^ + 2M^ - 3M^H). (30) 
To get 76 ^ 1, one needs to impose a less natural condition M^H"^ ^ M|, MfH. 

• 07 = 1/3 A?|7r(9,j7r)2 /a^ 
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where 

77 = {MlH'^) / {Ml^eH"^ + 2M4 - 3Mi3i7). (3I) 

The same consideration as for the case of M5 apply here. Note again that the numerical 
values, especially in the fifth and sixth benchmark points tend to be much larger for the 
M coefficients with the most spatial derivatives, thus conferming our expectations. 
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with 

78 = {M-jH^) / {Ml^eH'^ + 2M^ - 3MfH). (32) 

Notice that in this case , and in the following ones, the 7 coefficients are naturally ex- 
pected to be smaller than unity. For this interaction, as well as the two following ones, 
the numerical factor coming from the integration can be quite large, especially in the 
fifth and sixth benchmark points. This is clearly due to the six space derivatives that 
characterize these interactions. 

• 0, = -l/6 Msdfnid.kTiY /a6 
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with 



79 = {MsH^)/{Ml^eH'^ + 2M^ - 3MfH). 



(33) 



Let us stress that one could guess the amphtude of M^^q terms by simply looking at the 
results obtained for the interaction term tuned by Mj because, althought these terms 
might produce a different shape for non-Gaussianities, they have essentially the same 
structure as far as the integration is concerned. 
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-1/6 Mg dijTrdjkirdki'n- /c 
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where 



710 = {MgH^)/{M^ieH^ + 2M^ - 3M^H). 



(34) 



One could again just read off the maximum value from /^^, they differ by just a factor 
1/8. 



4-3. Some general considerations on the amplitudes of non-Gaussianities 

It might be worth now to pause and comment on our findings in relation to earlier results 
in the literature. As we have stressed, the non-Gaussianity generated from operators 
proportional to the M's masses arises from curvature terms that are often neglected. 
From general considerations, one expects that the curvature terms, especially the ones 
coming from SK^-tjpe of contributions, can often be neglected. This is not generically 
true though. The general structure of the amplitudes we have found for the bispectra 
can be schematically written as 



A 



1 




2M| + M^^eH^ - 3MfH 



(35) 



where M„ stands for a generic mass coefficient of the type M, M. 
If Mn 3> H/{ao + (ctg + 8/3o)^^^)^^^, and the masses are all of the same order, then it 
is natural to expect the M„ with the largest exponent to dominate. In this case the 
dominant terms are those associated with M2 and M3, corresponding to DBI inflation. 
Furthemore, we may also recall that large masses appear also in the definition of ao,/3o 
and, at least in DBI, a large M2 is needed to have a small speed of sound. 
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On the other hand, one can employ the freedom for all the M's not to be of the 
same order in magnitude. The M„'s have a natural upper bound that must be smaller 
than the Planck mass and in the general theory with both DBI and ghost (and more 
in general curvature terms) operators switched on, one can allow for a very small speed 
of sound without assuming much on the unconstrained masses. In such cases, the gen- 
eralized sound speed can be so small that H/{ao + {a^ + 8/3o)^^^)^^^ is actually larger 
than M. Consequently, the amplitudes of the masses with a smaller exponent, typically 
the barred M's, are not negligible any longer. We conclude that large non-Gaussianities 
may be induced by theories parametrized by suitable values of the M masses. 

It is important at this stage to offer some additional comments on the general 
structure of the Lagrangian and of the resulting amplitudes. As clear from Eq. ( l23l) . 
each space derivative acting on the scalar tt can be schematically written in Fourier 
space as a term Hn divided by a number which is much smaller than unity in the cases 
of interest. Now, it is clear from dimensional analysis that other terms in the action, 
specifically the ones with less derivatives, will have their own H/{aQ + {a^ + S/^o)'^^^)^''^ 
factors coming from derivatives but, most importantly, will also need some generic M 
factor to make the action properly dimensionless. This is exactly what happens in our 
theory (see Eq. ( !T3|l ) as well: the exponent of Mi is bigger than the one of M2 and 
so on. We have employed in the analysis of the amplitudes some freedom on these 
Mn parameters to show that even terms with higher spatial derivatives can give non 
negligible contributions to the amplitude. Of course, if one wants to have a reliable 
effective theory the M„ coefficients must eventually prevail over H/ (ao + (ao + 8/3o)"^''^)^^^ 
(in a Lorentz invariant theory it would suffice to ask for M„ ^ H, here we need more) 
so that it makes sense to consider higher derivatives up to some given finite order, but 
not further. 

5. The shapes of non-Gaussianities 

In this section we wish to analyze the shapes of the bispectra generated by the various 
operators analyzed above. In calculating the amplitudes on the bispectra in the previous 
section we have chosen to perform, albeit numerically, calculations with the exact 
wavef unctions. We decided to counterbalance the loss of information in not having 
aQ, Pq explicit in the result by running the same procedure for six different benchmark 
values of the parameters. In analaysing the shapes we do not enjoy the possibility 
to use the exact wavefunction any longer as we must approximate the wavefunctions 
inside the conformal time r integral(s) calculated as prescribed by the in-in formalism. 
Of course, exact results are always available when the classical solution reduces to the 
usual, Hankel function H{2,/2)ij) ■ In computing terms like 




(36) 
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we have choosen to expand in series each fc-mode inside the integrand within a region 
that starts from shghtly inside its effective horizon (setting to zero the function in the 
rest of the interval). This enables us to keep both, parameters and external momenta, 
arbitrary. This approximation is justified by the fact that, due to the oscillatory 
behaviour of the wavefunctions|| , the main contribution to the integral comes only 
from the region where all the wavefunctions are not oscillating anymore (to be safe, we 
actually choose to include the region where the /c-mode with the latest effective horizon 
is still within its horizon). Furthermore, it is reasonable to assume that, even if some 
non-negligible contribution is being left over because of this approximation procedure, 
it might have a systematic effect on all /c-modes and would therefore not change the 
shapes of the bispectrum of perturbations. This procedure has been calibrated with the 
calculations of Ref. [H] which were performed exactly: we are able to reproduce the 
very same shapes for the bispectrum. 

Following the scheme used for the amplitudes, we now report below the shape of 
non-Gaussianity for each interaction term. Within a single interaction term, we consider 
four configurations in the {a^, PQ)-p\a.ne. We preliminary found that the shape are not 
particularly sensible to their absolute values. This is expected as these coefficients, 
being the same for each fc-mode, enter mainly in the amplitudes, not in the shapes. 
Therefore, we limited our attention to the ao/fSo ratios. We probed the following 
cases: ao 7^ 0,/3o — (DBI); ao — 0, /3o 7^ (ghost); ao 7^ 7^ /3o (with the solution 
interpolating between DBI and ghost inflation and ao playing the dominating role in 
determining where is the effective horizon, we call this configuration A) and ao 7^ 7^ /3o 
(with the effect of (3q leading in the expression for the horizon, we call this configuration 
B). 

Let us briefly add some general preliminary comments. We need to point out 
that when we use the name DBI might generate some confusion. Indeed, we use it 
also to describe shapes due to interaction terms coming from curvature perturbations, 
such as Ml, 2,3....- What we mean here is that we can employ in the integrations the 
usual wavefunction Tik oc e~^^^^'^{l +^^/aokT) without resorting to any approximation. 
This wavefunction can be used as far as it is the solution to the equation of motion. 
This is a condition concerning only second order perturbations: we need only require 
Po = = Mo to employ it. We can, at the same time, have extrinsic curvature- 
driven interaction terms, the M's, and yet find exact results. The shapes corresponding 
to what we denominate DBI configuration have been thoroughly investigated in a 
number of papers [8], but what has been generally left out is the contribution coming 
from the extrinsic curvature terms: it is indeed possible to have the usual Hankel, 
-^3/2 (t), wavefunction as a solution to the classical equations and, at the same time, 

II Of course, an oscillatory behaviour is not, by itself, enough to provide a cancellation all over the 
-inside the horizon- region. Indeed in some cases for t — > —00 the amplitude of the wavefunctions 
increases and one certainly does not expect this to give zero contribution. On the other hand, much 
like in the simplest single field slow-roll calculations, one is expected to slightly rotate t to the imaginary 
plane to match the vacuum thus basically putting to zero the contribution at —00. 
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switch on curvature operators like M4,5,6. 7,8,9- All the shapes obtained in the DBI 
configuration are generated through exact analytical methods. We will call them 
as "exact-DBI configurations". Should we find, as we will, a shape which is not 
equilateral in the first configuration, that shape suffers none of any possible limitations 
the approximated method might introduce. The ghost configuration has been analyzed 
in depth in many articles, among which [331 ^S\, and again, not all the terms coming 
from extrinsic curvature have been taken into account. It is important to note though 
that in [231 133| 138] the M coefficients multiplying curvature terms have all been chosen 
of the same order thus resulting in only a couple of leading curvature (they correspond 
to the Ml, Mq contributions). Finally, the shapes in configuration A, B have never been 
analyzed before. 

The shapes reported below have been obtained by employing the shape function 
5(1, X2, xa), where Xi = ki/ki, which posseses the same /ci_2,3-dependence as the three 
point function. What is plotted exactly is xlxlB{l,X2,X3)/B{l,l,l) in the region 
satisfying X2 > xs > 1 — X2 [10]. Let us remind the reader that a shape is called 
local when it peaks for a small value of, say, X2, with xi = x^ = 1; equilateral when 
it peaks in the equilateral configuration xi = 1,X2 = 1 = X3 and is called fiat for 
squashed triangles with 1 = X2 +X3. In particular, we find, as detailed below, that some 
novel curvature-generated terms produce a fiat bispectrum which specifically peaks for 
X2 = X3 = 1/2. In presenting the shapes we follow the same order and organization we 
employed for the amplitudes: 

5.1. Shapes from DBI-like interactions and first two curvature- generated terms. 



• Oi = -2Ml7r{d,7rf/a^ 




Figure 1. DBI configuration on the left, obtained using exact methods; the 
approximated ghost shape on the right. 

Note that, although more or less sharply, all the four plots are peaked in the equi- 
lateral configuration. 
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Figure 2. A configuration on the left, B configuration on tfie right. 
• O2 = 4/3 M|7r3 




ij.O 6J(. 



Figure 3. exact DBI configuration on the left; approximated ghost shape on the right. 




Figure 4. A on the left, B configuration on the right for the Ms-driven interaction 
term. 

We obtain equilateral shapes in all four cases. It is somewhat expected that the 
general, interpolating solution employed in configurations A and B, will give qualita- 
tively the same plot, we have verified it it in these first two rounds of shapes. 
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. = -1/2 Ml{dinfd]n/a^ 




Figure 5. exact DBI configuration on tlic left; approximated gliost slrape on tlie riglrt. 




0.0 0.0 



Figure 6. A on tfie left, B configuration on the right for the Mi-driven interaction 
term. 

Also the bispectrum generated by this interaction term has an equilateral shape; 
the last two plots show, employing the general wavefunction, that also the interpolating 
models produce an equilateral shape. 

. = -1/2 M2/4 {bH{djTi){d,TTf + ATTd^djTidjTr) ja^ 

The various bispectra peak in the equilateral configuration. 
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1.0 




Figure 8. A on the left, B configuration on the right for the Mq interaction term. 



5.2. Shapes from curvature-generated novel interaction terms. 



• 05 = -2/3 MItt^^^tt /a^ 

With the term tuned by M4 we start including in our description the contributions 
that have so far been neglected in the literature. 




Figure 9. exact DBI configuration on the left; approximated ghost shape on the right. 
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Figure 10. A on the left, B configuration on the right for the M4-driven interaction 
term. 

Equilateral shapes for all configurations are obtained. 
• 06 = 1/3 Ml7r{df7iy /a^ 




Figure 11. exact DBI configuration on the left; approximated ghost shape on the 
right. 




0.0 k3 0.0 



Figure 12. A on the left, B configuration on the right for the Ms-driven interaction 
term. 

Again, all the four bispectra peak in the equilateral configuration. 



23 



• Oj = 1/3 M|7r(9,j-7r)2 /a^ 




Figure 13. exact DBI configuration on the left; approximated gliost sliape on tlie 
right. 




Figure 14. A on the left, B configuration on the right for the Mg-driven interaction 
term. 

This is one of the interesting novel curvature-generated terms that give rise to a flat 
shape (precisely the plot peaks at fci = 1, A;2 ~ 1/2 ~ k^)- Note that for a very similar 
interaction, namely the one generated by O^, we saw an equilateral plot. Here, deriva- 
tives combine to provide a different fc-dependent factor outside the integral. Writing 
in Fourier space the interaction term we obtain something proportional to k-i'kz^.i 
while here we obtain something like (/c2 ■ A^s)^. 
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= -1/6 Mridfnf /a^ 





Figure 15. exact DBI configuration on the left; approximated ghost shape on the 
right. 





Figure 16. A on tlie left, B configuration on the right for the M7-driven interaction 
term. 



All shapes peak in the equilateral configuration. 
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.0^ = -l/QNhd^n{d,knf M 




Figure 17. exact DBI configuration on the left; approximated ghost shape on the 
right. 




1.0 x3 1.0 



Figure 18. A on tlie left, B configuration on the right for the Mg-driven interaction 
term 

This is a second interaction term that produces, just as for O-j, a flat shape for the 
bispectra. Comparing it in Fourier space with our findings for Og, one can see that it 
is due to the way the spatial derivatives are combined. As shown in Figs. 19 and 20, 
also Oio gives rise to flat-shape bispectrum. We can see that the interactions Cs, Oq 
and OiQ have the same structure as far as the integral is concerned; on the other hand 
their /c-dependence goes like kl k^, /cf (/c2 • k^Y + perm and (fci ■ A;2)(/c2 ■ k^){k^ ■ ki), 
respectively. The last two produce a flat shape. 
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Figure 20. A on the left, B configuration on the right for the Mg-driven interaction 
term. 

The shapes are peaked in the flat configuration. 

Obviously, having presented all the shapes due to each individual interaction, one 
might proceed with the study of the shape of linear combinations of them, much in 
the spirit of the orhtogonal shape recently introduced On the other hand, we 

are using approximated methods for three of the four configurations under scrutiny 
and it is therefore not a cautious step to infer new shapes from linear combinations of 
approximated ones, especially when delicate substractions are involved. 

5.3. Some general considerations on the shapes of non-Gaussianities 

From the above discussion we can read off some general qualitative features 
characterizing the shapes of the bispectrum of a generic single-field model of infiation: 

• Two qualitatively very different kinds of shapes appear: equilateral and flat. As 
pointed out before, the next natural step would be to consider the shapes resulting 
from linear combination of the various interaction terms in the spirit of Ref. [38] . 

• In single fleld-models of inflation usually large non-Gaussianities are associated 
to equilateral shapes. In order to obtain a flat shape one needs to consider linear 
combinations of interaction operators such as what is done in [38] or models with an 
initial vacuum different from Bunch-Davies. Interestingly, in our case a flat shape 
emerges simply from individual operators generated by curvature-related terms. 
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When the flat shape appears, it does so in all four conflgurations considered. We 
stress this point because it implies that this result does not depend on the type 
of wavefunction one employs in the calculation, be it the DBI inflation solution, 
the ghost inflationary one and in the intermediate cases. The results of the DBI 
conflguration are exact and easily reproducible with analytical methods. In fact, 
the DBI wavefunction is the usual solution of the standard single-fleld slow-roll 
inflation. Adopting such a wavefunction we may provide analytic results for the 
flat bispectra 

_{60a%'^ + lla%^ + b^ + {48a^-4a%^-3b^)c^-4b^c^)H^Mi 

(15a3 + 63 + 35c2) (I2a6 + fe^ _ Qb^c"^ + gb'^c^ + 8a^ {b^ - 2bc^)) H^M^ 



{CklCk2Ck3) 



3 (15a3 + 63 + 35c2) (Sa^ + 6^ - 6&^c^ + 8fe^ c^ + 8a3 (63 _ 2bc')) H^Mq 

2a? al 66 M^^ 



(CfclCfc2Cfc3)M9 = 0.9^,3^6.3 />^6 ' 

(37) 

where the overall momentum conservation delta has been omitted, 

a = {kik2hf'^] b = ki + k2 + h] c= (/ci/cs + /ci/ca + k2k:if/^, 

and ao is now the usual speed of sound: 

ao = -Ml^H/{-Ml^H + 2M^) = . 

Note that all the three expression given above have a maximum precisely in the flat 
conflguration {ki = 1, k2 = 1/2 = ks). What immediately stands out in Eq. fl37j) is 
the presence in the numerator of factors consisting of subtractions between generally 
positive fc-symmetrized terms: it is this characteristic that selects a flat, rather 
than a equilteral shape, as one can readily verify by checking the bispectrum of the 
"equilateral" interaction terms. 

The expression in Eq. fl371) is exact in the DBI case and can be employed to get a 
shape qualitatively similar in the other three conflgurations. For practical purposes, 
we give below a very simple expression that very closely mimics the behaviour of 
the typical bispectrum contribution that generates flat shape 

, , , {-kl + - k,^) {kl + k2' - k^) {-kl + k2' + k^) 
{CklCk2Ck3)M fc3fc2^fc33 {k, + k2 + k-s) 6 

As a general rule, the terms which are going to generate a flat shape can be read 
off already at the Lagrangian level: indeed the flatness originates from the way 
the external momenta combine with each other and are summed over. Whenever 
mixed space derivatives act on a single vr term and the mixing is repeated on at 
least another tt fleld, the shape turns out to be flat (note that this criterium puts 
-^6,8,9 contributions in the same, "flat" class, but correctly excludes apparently very 
similar ones such as M2,3,7)- 
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6. Conclusions 



The purpose of this paper was to employ a powerful tool such as effective field theory to 
obtain as general as possible a knowledge of primordial non-Gaussianities generated in 
a very general set-up of single field models of inflation. We have performed a study of 
the corresponding bispectra, in terms of both amplitudes and shapes. We extended the 
results existing in the literature in different ways. First, we have improved the treatment 
of the wavefunction describing the behaviour of the cosmological perturbations at second 
order in perturbation theory. Secondly, we computed the amplitude and the shape of 
bispectra for theories which interpolate between the most studied ones, i.e., the DBI and 
the ghost models. Third, we have pointed out the importance of the curvature-induced 
operators. Their study has revealed that large non-Gaussianities may be generated 
with a fiat shape which is quite uncommon for single-field models of inflation so far 
analyzed. It would be interesting to identify in which (class of) theories of inflation 
such operators arise. A natural way to extend our findings is to study with the same 
philosophy the four-point correlator, the trispectrum, and to identify, for instance, those 
classes of inflationary models where the bispectra are suppressed, but large trispectra 
are generated. These issues are currently under investigation [32]. 
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